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Abstract 



Four-dimensional Einstein-Maxwell-dilaton-axion system restricted to space- 
times with one non-null Killing symmetry is formulated as the three-dimensional 
gravity coupled sigma-model. Several alternative representations are discussed 
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isometries on the initial set of (non-dualized ) variables is found. New mulicenter 
solutions are obtained via generating technique based on the formulation in terms 
of the non-dualized variables. 
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1 Introduction 



Einstein equations in General Relativity have a nice property to admit a three-domensional 
gravity coupled sigma-model representation, being restricted to space-times possessing a 
non-null Killing symmetry. Moreoer, it turns out that the target space of this sigma-model 
is a symmetric Riemannian space, namely, the coset SL(2,R)/SO(2). This sigma-model 
may serve a convenient starting point for a number of solution-generating methods, most 
notable being an inverse scattering transform and Backlund transformations, which arise 
when further Killing symmetry commuting with the first one is imposed [[![]. The existing 
knowledge of a great variety of analytic solutions to the Einstein equations is mostly due to 
this property 0. 

Similar property is shared by other theories related to General Relativity: Einstein- 
Maxwell (EM) theory ||, higher dimensional vacuum Einstein equations [ffl, bosonic sectors 
of some dimensionally reduced supergravity theories ||, |J, @. Crucial feature of all these 
theories is the symmetric space property of the corresponding sigma-model in three dimen- 
sions (for a review see, e.g., Breihtenlohner and Maison ||). The symmetric space property 
in its turn requires that the target space possess a sufficient nubmer of isometries. This 
convenient way to investigate the problem is due to the idea of potential space introduced 
by Neugebauer and Kramer ||. 

Recent iterest to such models is motivated by the string theory which is likely to provide 
a consistent description of gravity at the Planck scale. The most promising string model con- 
sidered so far is that of the heterotic string. In the low-energy (classical) limit of this theory 
with compactified extra dimensions one gets in the bosonic sector the Einstein gravity cou- 
pled to massless vector and scalar fields. The simplest model of this kind — "dilaton-axion 
gravity" - - incorporates basic features of the full effective action, and can be formulated 
directly in D = 4 where it includes one U(l) vector and two scalar fields coupled in such 
a way that the theory possesses non-abelian SL(2, R) duality symmetry. This theory was 
a subject of recent investigations |§ which are summarized and developed further xin the 
present paper. 

The Einstein-Maxwell-Dilaton-Axion (EMDA) coupled system contains a metric g^, a 
U(l) vector field A^, a Kalb-Ramond antisymmetric tensor field B^, and a dilaton (p: 



fj,, is, A = 0,1,2,3, and oj^l, ^3YM are Lorentz and Yang-Mills Chern-Simons three-forms 
respectively. In four dimensions a Kalb-Ramond field is equivalent to the Peccei-Quinn 
pseudoscalar axion k, 



16vri I 




where 



H^x = d^B vX - A^F vX + cyclic, 



(1.2) 



" 2 dx T ' 

and the corresponding Legendre transformation of the action (1.1) reads 




(1.3) 



s = wJ{- R+ 




where = \E^ Xt F Xt . 



Our purpose is to investigate hidden symmetries associated with this model when an 
additional assumption of stationarity is imposed on the space-time metric. The plan of the 
papers is as follows. In Sec. 2 we outline the derivation of a three-dimensional sigma-model 
from the stationary EMDA system. In Sec. 3. complex potentials are introduced which make 
formulation similar to the Ernst formulation of vacuum and electrovacuum theories. The 
corresponding matrix formulation is given in Sec. 4. An alternative description in terms of 
non-dualized quantities is then suggested (Sec. 5) and used to genetrate EMDA multicenter 
solutions (Sec. 6). 

2 a— model in three dimensions 

The system of equations corresponding to (1.4) consists of the dilaton-axion modified Maxwell 
equations 

Vvie'^F^ + kF> m/ ) =0, (2.1) 



dilaton and axion equations 



V M VV = l -e~ 2 *F 2 + l -e A \dK)\ (2.2) 



V M (e 4 V"U«) + = 0, (2.3) 

and the Einstein equations 

iV = 2^,„ + \e^^K v + e- 2 ^(2F lxX F\ + ^F 2 g fiu ). (2.4) 

Here F 2 = F^F^", (dn) 2 = g^ il, {d IJi K){d v K), and is the covariant derivative with respect 
to the 4-dimensional metric g^ u . 

Consider a space—time admitting (at least) one Killing vector field which we choose to 
be timelike. Then it is standard to write an interval using the three-metric hy, the rotation 
vector u>i, = 1, 2, 3) and the scalar /, depending on the space coordinates x\ as follows 

ds 2 = g^dx^dx" = f(dt - oOidx 1 ) 2 - jh ij dx t dx j . (2.5) 

To derive a three-dimensional a-model one has to introduce the set of appropriate vari- 
ables. It consists of the electric v and magnetic u potentials as well as the twist potential 
related to u>i through a certain curl equation. Electric and magnetic potentials are introduced 
through the relations 

F t0 = -j=diV, (2.6) 

e -2<t, F ij + K pij = f e ijk dkU (2.7) 
y2h 



The first substitution solves the spatial part of the Bianchi identity 



V M F^ = 0, 



(2.8) 



while the second satisfies the spatial part of the modified Maxwell equations (2.1). From 
(2.6) and (2.7) the remaining Maxwell tensor components can be expressed through v, u, f, uJi 
and hij. Introducing instead of u>i a three vector 



e ijk 

T l = -f 2 -j=djU k , 



(2.9) 



one obtains from the remaining Maxwell equations the following set of three-covariant equa- 
tions 



Vh f P Vh f 



je 2 *Vhh ij (d jU - Kd jV ) 



- ^t%v = 0. 
f 2 



0, (2.10) 



(2.11) 



The equations for dilaton and axion fields (2.2), (2.3) with account for (2.5) will take the 
form 



1 



1 



^di(Vhti j d j( j)) = — e'^diV&v - e^fim - ndivX&u - nd l v) 
Vh 2 j 

-^diiVhti'e^diK) = -e 2<t> (diu - nd^&v, 
Vh J 



e 4</> 

+ —diKd'n, (2.12) 



(2.13) 



where (and in what follows) raising and lowering of the three-indices is understood with 
respect to the three-metric and its inverse h l K 

Remarkably, the mixed components of the Einstein equations (2.4) remain unaffected by 
neither a dilaton nor an axion: 



i 



iVh 



e ijk 
'h 



e tjk d k Tj = —j=djvd k a. 



This equation can be integrated giving the twist potential x 

n = diX + vdiU - udiV. 



(2.14) 



(2.15) 



An equation for x can De found by taking a three-covariant divergence of (2.15) multiplied 
by f~ 2 with account for the definition (2.9): 



where 



/(A% + vAu — uAv) = 2(diX + vdiU — udiv)d l f, 



A = J-di(y/hti j dj) 
Vh 



(2.16) 



(2.17) 



is a three-dimensional Laplacian. Finally, an equation for / follows from the (OO)-component 
of the Einstein equations in a similar form 



fAf - dtf&f + nr 1 = f le^idiU - Kdiv)(&u - nd l v) + e-^divtf 



(2.18) 



where the definition (2.15) has to be used. 

The system of 6 equations (2.10)-(2.13), (2.16), (2.18) constitutes the set of equations 
for a curved space a-model with the 6-dimensional target space. The remaining (spatial) 



Einstein equations from (2.4) may be consistently regarded as 3-dimensional Einstein equa- 
tions for the metric hij. Their source part is derivable from the same sigma-model as the 
3-dimensional energy momentum tensor 

Kij = Jpififj + T iTj) + 24> il (j) j + -e H KiKj- 

- j e'^v-jVj + e 2<t, (u ti - Kv,i)(u tj - kvj) , (2.19) 

where r, is given by (2.15), and IZij is the three-dimensional Ricci tensor. 

Now, it can be easily checked that all 3-dimensional equations may be obtained by a 
variation of the action 

S = J{U- ^[(V/) 2 + (Vx + vVu - uVvf] - 2(V$) 2 - 

- \e^(VK) 2 + i [e 2<t> (Vu - kVv) 2 + e^^vf] }Vhd 3 x, (2.20) 
2 / 

where 1Z = 1Z\ and V stands for 3-dimensional covariant derivative. 

This action can be rewritten as the gravity coupled three-dimensional a-model 

S = J (n - g A BdiV A d 3 tp B h ij ) y/hd 3 x, (2.21) 
where <p A = (/, Xi v i u i K i 4>i) A — 1> ■■•,6. The corresponding target space metric reads 
dl 2 = -r 2 [df + (d X + vdu - udv) 2 ] - f-^idu- ndvf + e~ 2(t, dv 2 } 

+ 2d(j) 2 + ]^e H dK 2 . (2.22) 

This expression generalizes Neugebauer and Kramer potential space metric found for the 
stationary Einstein-Maxwell system in 1969 0]. It is worth noting that the present a-model 
does not reduce to the Einstein-Maxwell one if k — (j) — 0, because the equations for a 
dilaton and an axion generate constraints F 2 = FF = 0. Hence, generally the solutions of 
the Einstein-Maxwell theory with one Killing symmetry are not related to solutions of the 
present theory by target space isometries (except for the case F 2 = FF = 0). From the other 
hand, one can consistently set in the present a-model / = x = = K = O reducing it to 
the Einstein vacuum sigma-model. Therefore all solutions to the vacuum Einstein equations 
with one non-null Killing symmetry are related to some solutions of the system in question 
by target space isometries. 

3 Complex potentials 

More concise formulation of the theory can be achieved in therms of the following complex 
variables 

z = K + i e - 2 <t> (3.1) 
$ = u - zv (3.2) 
E = if- x + v$ (3.3) 



The first is the standard complex dilaton-axion field, while two other may be regarded as 
suitable generalizations of the Ernst potentials. In terms of them the target space metric 
reads 



dl 2 



1 

2P 



,„ 2Im$ ^ 

dE + rf$ 

lmz 



/Im$ V 



V Imz / 



dz 



f\mz 



Im$ 

\mz 



dz 



+ 



\dz\ 



2(lmz) 



(3.4) 



It can be shown that this three-dimensional complex manifold is Kahler. Indeed, in terms 
of the complex coordinatess z a = (E, —z, $), z a = (E, —z, $), where z a is a shorthand for 
z a , one can write an Hermitian metric on the target space as 



dl 2 = K a - p dz a dzP. 
This implies the existence of a non-degenerate two-form 

n = l -K a - p dz a A dz?. 

One can easily check that this form is closed: 

dfl = 0, 



(3.5) 



(3.6) 



(3.7) 



what means that we have a Kahler manifold. 

Two describe isometries of the target space it is convenient to take into account the 
following two dicrete symmetries. The first is quite obvious 

z = E, E' = z, $' = $. (3.8) 

Another useful discrete symmetry operation, which will be called "double prime", is 

Z n E -/ $ 



E 



Ez + <$> 2 



z 



Ez + <$> 2 ' 



Ez + $ 2 ' 



(3.9) 



The set of coninuous isometries consists of three gauge (electric e, magnetic m, and grav- 
itational g) transformations, a scale transformation (s), SL(2,R) duality (^1,^2,^3), two 
Harrison-type (Hi,H 2 ) and Ehlers-type (E) transformations. The Killing vectors corre- 
sponding to gauge transformations in terms of complex potentials read 



K e = 2$d E - zdq> + c.c, 



■ <9$ + c.c, 
d E + c.c. 



The scale generator is 



K s = 2Ed E + $<9$ + c.c. 



(3.10) 

(3.11) 
(3.12) 

(3.13) 



The Ehlers-Harrison-type generators can be obtained by "priming" the gauge gener- 
ators 

K Hl =K' e = 2§d z - Ed* + c.c, (3.14) 
(Ez — $ 2 )<9$ — 2§(zd z + Ed E ) + c.c, 



K 



H 2 



k: 



K e = K" = $ 2 c), - E(Ed E + $d*) + c.c. 



(3.15) 
(3.16) 



Two generators of the ^-duality subgroup correspond to the primed g and s transforma- 
tions 

K dl =K' g = d z + c.c, (3.17) 
K d3 = K' s = 2zd z + + c.c, (3.18) 
while the remaining one may be obtained by priming the Ehlers generator 

K d2 =K E = § 2 d E - z{zd z + $c\) + c.c. (3.19) 

It can be checked that these 10 Killing vectors form the algebra sp(4, R), while the target 
space in isomorphic to the coset Sp(A, R)/U(2). Finite Sp(A,R) transformations may also 
be obtained starting with the gravitational 

E = E + X, $ = $ , z = z , (3.20) 

electric 

E = E - 2A$ - \ 2 z , $ = $o + Az , z = z 0: (3.21) 

the magnetic 

E = E , $ = $ + A, z = z , (3.22) 
gauge as well as the scale transformation 

E = e 2X E , $ = e A $ , z = z , (3.23) 

(which can be easily integrated from the infinitesimal form). The remaining six group ele- 
ments can be found using prime operations. The 57(2, R) -duality subgroup read 

= $o, z = z + A, 
+ \z 1 + \z 1 + \z 





E 


= E , $ 


d 2 : 


E 


= E + A- 


d 3 : 


E 


= E , $ 



s A $ , z = e 2X z , 
the electric and magnetic Harrison transformations are 

E = E , $ = $ + \E , z = zq- 2A$ - X 2 E , (3.25) 
E = Eo ^ _ %(1 + \<5> ) + \E z z 



(1 + A$ ) 2 + \ 2 E z ' (1 + A$ ) 2 + \ 2 E z ' (1 + A$ ) 2 + \ 2 E z ' 

(3.26) 

while the Ehlers-type transformation reads 

E ~TTW^ $ -TTA£o' z - Zo + TTW, (3 ' 27) 

The isometry group acts transitively on the target space which is a symmetric Rieman- 
nian space. The direct way to show this is to introduce matrix representation in terms of 
symmetric symplectic 4x4 matrices forming a coset Sp(A, R)/U{2). 



4 Matrix representation 

For any matrix G G Sp(A, R) one can perform a Gauss decomposition 

G = G L G S G R , (4.1) 

where 

- (V '.')• c(i m 

and S, R, L are real 2x2 matrices, R, L being symmetric, R T = R, L T = L. Useful relations 
are 

Gr 1 Gr 2 = Gr 1+ r 2 , Gs 1 Gs 2 = Gs 1 s 2 , Gl x Gl 2 = Gl 1 +l 2 , (4-3) 

so that for inverse matrices one has G^ 1 = G-r, Gg 1 = S-r, G^} = G-l- With this 
parametrization 

S T ~ l S T ^R ' 

LS 7 " 1 S + LS^R 



G=[ r Vi Cl " I- (4.4) 



A subset of symmetric symplectic matrices M e Sp(4,R), M T = M represents a coset 
Sp(A,R)/U(2). The corresponding Gauss decomposition reads 

M = Q T VQ, (4.5) 

Q-(S«). M 

where two real symmetric 2x2 matrices are introduced Q T = Q, P T = P. Multiplying 
matrices we obtain 

M =(<P p + qEq)- < 47 > 

It is useful to combine P and Q into one complex symmetric matrix 

Z = Q + iP, (4.8) 
which transforms under an action of Gr, Gl, Gs as follows: 

Z ^ Z + R, (4.9) 

for G = G L : 

Z'^Z^-L, (4.10) 

while for G = G s ■ 

Z^S T ZS. (4.11) 
Here the following transformation law for M is assumed 

M -> G T MG. (4.12) 

These transformations are precisely the matrix analog of the SL(2, R) transformations 
which hold for the usual scalar dilaton-axion field (3.1) under S-duality generated by 
Kd 1} K d2 , K ds . 



It can be shown that the explicit expression for Z in terms of the complex potentials 
introduced in the previous section is very simple 

Mf-!)- (413) 

Sigma-model of the Sec. 3 can now be rewritten as a coset model for the matrix M (4.7). 
Let us introduce a matrix current 

J M VMM '. (4.14) 

Using (4.8) into we get 

M= ( -P^VP-P^VQP^Q P^VQP- 1 ^ 

~\VQ- QP- l VP - VPP' l Q - QP^VQP^Q VPP' 1 + QP^VQP- 1 J ' 

In this notation, the three-dimensional sigma-model action can now be written as 

S = J (-K + l -TriJ M2 )^j Vhd 3 x, (4.15) 

while the corresponding field equations read 

VJ M = 0, (4.16) 

n mn = -Tr{J M rn J M n ). (4.17) 

There exists also a concise representation of the same model directly by 2 x 2 matrices. 
Indeed, from three 2x2 matrix equations following the present representation 

V (P^VQP- 1 ) = 0, (4.18) 

V (QP^VQP- 1 + VPP' 1 ) = 0, (4.19) 

V (VQ - VPP- l Q - QP- l VP - QP^VQP^Q) = 0, (4.20) 
only first two are independent. The last can be rewritten using first two as 

PV (P^VQP" 1 ) P = 0. (4.21) 

Introducing two 2x2 matrix currents 

J x = VPP' 1 , J 2 = VQP- 1 . (4.22) 

we can rewrite the equations of motion as 

VJi + J 2 2 = 0, 

VJ 2 - Ji J 2 = 0, (4.23) 

T^mn — {JlmJln + </2m^2n) • 



It can also be checked that 

Tr (j x 2 + Jf) = ^TrJ 2 . (4.24) 

Further simplification comes from the introduction of a complex matrix dilaton-axion Z 
instead of its real and imaginary parts (from an explicit for of Z ) it is clear that real and 
imaginary parts are more involved because of the complicated structure of the denominator). 
Then instead of two real 2x2 currents one can build one complex matrix current 

J z = VZ(Z-Zy\ (4.25) 

and cast the action into the following form 

S = (-TZ + 2TrJ z J^) Vhd 3 x. (4.26) 
The correspoding equations read 

VJ Z = J Z (J Z + J Z ), 

K mn = 2TrJ Zm J Zn . (4.27) 

The metric of the target space can be presented in terms of the matrix-valued one-forms 
corresponding to the above currents 

to = dMM~\ u! = dPP-\ u 2 = dQP-\ u z = dZ(Z - Z)~\ (4.28) 

by three alternative ways 

ds 2 = ]Tr{u 2 ) = ^Triuj 2 + uj 2 2 ) = 2Tr{u z u z ). (4.29) 

5 Back to non-dualized variables 

To achieve sigma-model description of the system in three dimensions it was essential to use 
the dualized variables for the rotational part of a metric and for the magnetic part of a vector 
field, as well as to introduce the Peccei-Quinn pseudoscalar counterpart to the antisymmetric 
Kalb-Ramond field. This description was particularly useful in deriving hidden symmetries 
of the dimensionally reduced system. Once they are found, one can go back to the initial 
formulation and to reveal their action on the non-dualized set of variables. Apart from 
purely formal interest, this also opens a way to obtain an infinite symmetry algebra in the 
further two-dimensional reduction. 

Transition to the initial variables may now be regarded as the formal solution of the 
matrix equation 

VJ M = 0, J M = VMM- 1 . (5.1) 

Since its divergence is zero, the matrix current can be presented as a curl of some vector 
matrix iV 

V x N = J M . (5.2) 
It can be checked that an alternative form of the equation VJ M = reads 



(5.3) 



In terms of 2 x 2 symmetric matrices P ( Q one gets the following set of equations 



V x Q = P^VQP -1 , (5.4) 

V x U = VPP- 1 + QP- 1 VQP-\ (5.5) 

V x V = VQ - VPP- l Q - QP- l VP - QP^VQP^Q, (5.6) 
where Q, V and U are symmetric real 2x2 matrices forming N as follows 

From explicit expressions for P and Q 

p _ / / - e^v 2 -e'^v \ 



wv-x w 

y w —k j 

it is clear that they provide a natural decomposition of the full set of sigma-model coordinates 
into non-dualized (/, 0, v) and dualized (%, w, re) subsets. Hence, it can be anticipated that 

— * 

the matrix P will remain unchanged, while Q will be replaced by Q. After some calculations 
one can express it through the non-dualized variables as follows: 

Q-( Q -(A + A u;) \ 

\-(A + A uj) -B + A (A + A u;) )' ^ 

where A = v,P>i = B 0i , ui is the metric function and Ai is three-dimensional vector- 
potential defined as \f2F^ = (diAj — djAi). Two other matrices U and V can be found 
once Q is known (all three-dimensional vector operations in the above formulas have to be 
effected using the metric h lJ . 

It is worth noting that the spatial components of the Kalb-Ramond field do not enter 
into the above set of variables. In fact, in the stationary case they can be expressed through 
other variables. Introducing a three-dimensional dual to 

C l = ^ k B jk , (5.9) 
from stationarity re = one can find after some rearrangement 

VC = V(B xu) + [B- A (A + A lo)]V x Q + (A + A Q u)V x (A + A O). (5.10) 

One can derive the three-dimensional action variation of which gives the equations for 
matrices P and Q. First, from the Eq.(5.4) one gets 

V x [P(V x n)P] = 0, (5.11) 

while (5.5) takes the form 

V(VPP _1 ) + P(V x Q)PV x n = 0. (5.12) 



Together with the Einstein equations 

R l3 =2Tr(jfjf + J?J?), (5.13) 

where 

J p = VPP -1 , J n = PV xtt (5.14) 

we have the complete description of the system in terms of P and Q. The Eqs. (5.11-12) 
can be rewritten as 

VJ P + (P 1 ) 2 = 0, V x J n - J n x J p = 0, (5.15) 
while the corresponding action read 

S = J d z xh^(-R + 2Tr[(J p ) 2 - (J n ) 2 }). (5.16) 

The equations for U and V still contain undesired matrix Q. However two "Legendre 
transformed" quantities 

U' = U-Qtt, V' = V + QVLQ + U'Q + QU' t (5.17) 

— * 

satisfy the equations depending only on P ( Q: 

v x u' = vpp- 1 - p(v x n)p x (i (5.18) 

V x V' = P(V x Q)P + P(V x Q)P x f7' T - [7' x P(VQ)P. (5.19) 
Similarly, a non-dynamical variable C may be replaced by 

C' = c- \[B + A {A + A uj)] x u, (5.20) 

the equation for which contains only "allowed" variables: 

VC' = iTr(fieV x fie). (5.21) 

where 

-(-".:) 

— * 

Target space isometries induce certain transformations of the matrix f2. These can be 
found by a direct computation using the Eqs. (5.4-6). The subgroup R leaves VP and 
VxO unchanged, and now some more general transformation of the same kind can be foun 

P^P, tt^tt + VR', (5.22) 

where R' is an arbitrary non-constant symmetric matrix. 
Now, for the ^-transformation it is convenient to redefine 



-IT 



(5.23) 



so that we get the following counterpart 

P- 1 -> S^P^S', (5.24) 

n -> S' T QS'. (5.25) 
Finally, for the L-transformations one has 

P" 1 = P _1 + LQoPo 1 + ^cT'Qo^ + £(P + QoPo'Qo)^ (5.26) 

while f2 can be shown to transform as 

(l -> + LC7 + f7 T L + LT?L. (5.27) 

For generating purposes it is convenient to combine the action of S, L transformations. 
Omitting primes, one gets 

P- 1 = S T [P 1 + LQoP' 1 + P^QqL + L(P + Q P^Q )L]S (5.28) 

Q — S T [Qo + LU + UqL — LV L]S. (5.29) 

It is worth noting that for the three-dimensional description in terms of the non-dualized 
quantities we have introduced more variables than in the previous sigma-model formulation. 
The excessive variables could in principle be eliminated by fixing the gauge, but then three- 
dimensional covariance would be lost. 



6 Multicenter solutions 

Here we give some examples of the application of symmetries to derive multicenter solutions 
to dilaton-axion gravity. Le us specify the transformation formulas for asymptotically flat 
solutions with zero asymptotic values of material fields. In this case 

P o£> = a 3, Qooc = ^Ooo = 0. (6.1) 

Demanding that the L-transformation leave unchanged asymptotic values of /, 0, v, one 
obtains for the matrix L the following expression 

L = IE, (6.2) 

where 

and o — il, I — const. We will assume that the S'-transformation leads to a non-trivial 
asymptotic for the dilaton 

^ -e 2 *°° J ' 
These assumptions leave one free parameter in S: 

/ ch6 e*°°sh6 \ 
[she e*°°ck9 j ' 



As a whole we get two parameter (/ and 6) family of solutions applying S, L-transformation 
to a chosen asumptotically flat seed solution. 
As a first example we take as a seed 

Pq- 1 = <7 3 , Q = AS' (6.3) 

— * 

where E' corresponds to the replacement of a to a' in (). A (non-constant) matrix vector A 
is introduced according to 

N \ 

VxA = VA, A = ]T— (6.4) 

where r*j are positions of the sources. Then three-metric turns out to be flat 

dl 2 = dr 2 . (6.5) 

Physically this solution corresponds to an equilibrium system of an arbitrary number of 
point-like magnetic monopoles endowed with axion charges (magnetic repulsion is balanced 
by axion attraction). This field configuration generates twist, so that the resuling metric is 
a massless multi-Taub-NUT. 

Now we apply to this solution a combined S, L transformation described above () setting 
a = —a (only in this case L-transformation gives non-trivial result). One gets 

p- 1 = (i + 4ja) ( i 

u ~ A y e {-ae*- e 2 ^ ) ^ 6 { ae*" e 2 ^ ))• 
In physical terms 

f- 1 = 1 + 4ZA, e 2<t> = e 2 ^f-\ A = (6.6) 
Q = ( e - 2ae - Al 2 e 2ae )K A = ae^ {e' 2ud + Al 2 e 2ae )A, B = -e 2 ^u. (6.7) 
In the case 

after redefinitions 21 X — > A, crl/\l\ — > a the solution can be written as 

f-i = e 2(^oc) = i + 2X, A = 2ae^A, A = u = B = 0. (6.9) 

and therefore describes the system of point massive magnetic charges endowed with dilaton 
charges. In this case (no axion field) twist is not generated. 

As a second example consider an (S, L)-transformation of the electric Majumdar-Papapetrou- 
like solution 

P = (7 3 + AE, tt = 0, dl 2 = dr 2 , (6.10) 

corresponding to 

Q = n = V = 0, (6.11) 



With no additional assumptions for a ( a, one obtains 



0-9/A (I 1 ^ ' I ch2d e* a °sh26 \\ 
" ~ ^ -e 2 ^ J + °° \ e+°°sh26 e 2 ^ch29 ) ) ' 

Extracting physical parameters, one can see that we have generated magnetic, axion and 
NUT charges. 
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